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1 Introduction 



Integrable models in two-dimensional quantum field theories and the integrable vertex models in two 
dimensional classical statistical mechanics are the common link of many relevant exactly-solved quantum 
models in one dimension [1, 2, 3, 4]. Examples are the XXX, XXZ and the XYZ Heisenberg chains that 
find more and more applications in contemporary physics. The quantum inverse scattering method or 
algebraic Bethe ansatz is the tool to exploit this fact systematically [3]. The key toward this synthesis 
is the existence of a central object 7^(w), where u is a spectral parameter, acting on the tensor product 
y (g) y of a given vector space V and being a solution of the celebrated quantum Yang-Baxter equation 

TZi2{u)TZi3{u + v)n2s{v) = 7^23(^;)7^l3(w + ^;)7^l2(u), (1.1) 

in yi y^ y^, where ni2='R-(g) 1, 7^23 = 1 «) 7^, etc. 

The solution 7l{u) of (1.1) is said to be quasi-classical if TZ{u) also depends on an additional parameter 
T] in such a way that 

TZ{u,r]) = l+r] r{u)+o{r]^), (1.2) 
where 1 is the identity operator on the space V ^V. The "classical r-matrix" obeys the equation 

[ri2(u),ri3 (« + !') +r23{v)] + [nsiu + v) , r23{v)] = 0. (1.3) 

This equation, called the classical Yang-Baxter equation, plays an important role in the theory of classical 
completely integrable systems [5]. 

Nondegenerate solutions of (1.3) in the tensor product of two copies of a simple Lie algebra g , 

rij{u) e gi (Xi gj , i,j = 1, 2, 3, were classified by Belavin and Drinfeld [6]. 

Gaudin models [7] constitute an important class of one dimensional many-body systems with long- 
range interactions with important applications in many branches of fields ranging from condensed matter 
to high energy physics [8, 9, 10, 11, 12]. They are related to the classical r-matrices through the definition 
of the density Gaudin Hamiltonians 

JV 

Ga^'^VabiZa- Zb) (1.4) 
bj^a 

The commutativity condition [Ga, Gb] = is a consequence of the classical Yang-Baxter equation (1.3). 

The classical Yang-Baxter equation also has a counterpart in conformal field theory, which can be 
described in the following way: in the skew-symmetric case rj,(— u) -|- rij{u) = 0, it is the compatibility 
condition for the system of linear differential equations 

d'i/{zi,...,ZN) 



dzi 



'^nj{Zi- Zj)^{zi,...,ZN) (1.5) 
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in N complex variables zi,...,zn for vector- valued functions 4* with values in the tensor space V = 
(g) • ■ • , where k is a coupling constant. 

In the rational case [6], very simple skew-symmetric solutions are known: r{u) = C2/U, where C2 € 
g ® g is a symmetric invariant tensor of a finite dimensional Lie algebra g. The above system of linear 
differential equations (1.5) is known as the Knizhnik-Zamolodchikov (KZ) system of equations for the 
conformal blocks of the Wess-Zumino-Novikov-Witten (WZNW) models of the conformal field theory on 
the sphere [13]. 

The work of Babujian and Flume [14] unveils one link between the algebraic Bethe ansatz for the 
theory of the Gaudin models and the conformal field theory of WZWN models. In their approach, the 
Bethe wave vectors for an inhomogeneous lattice model give, in the quasi-classical limit, solutions of the 
KZ equation for the case of simple Lie algebras. For instance, in the su(2) example, the quantum inverse 
scattering method [15] allows one to write the following equation 

t{u\z)^{ui^... Up) = A{u, Ui, - ■■ , Up\z)^{ui, • • • , Mp) - • (1-6) 

a=l 

Here t{u\z) denotes the transfer matrix of the rational vertex model in an inhomogeneous lattice. 
meaning = $(«!,••• u, Ma+i, Wp) ; ^a('Ui,--- ,Up\z) and A(u, ui,-- - ,Up\z) being complex 
valued functions. The vanishing of the so-called unwanted terms, J^a = 0, is demanded in the usual 
procedure of the algebraic Bethe ansatz [16, 17] in order to determine the eigenvalues of t{u\z). In 
this case the wave vector $(ui, • • • ,Up) becomes an eigenvector of the transfer matrix with eigenvalue 
K{u, ui, • • • ,Up\z). If we keep all unwanted terms, i.e. J^a 7^ 0, then the wave vector $ in general satisfies 
the equation (1.6), named in [18] as off-shell Bethe ansatz equation (OSBAE). 

There is a close relationship between the wave vector satisfying the OSBAE (1.6) and the vector-valued 
solutions of the KZ equation (1.5): the general vector valued solution of the KZ equation for an arbitrary 
simple Lie algebra was found by Schechtman and Varchenko [19]. It can be represented as a multiple 
contour integral 



^{zi, . . . , ^jv) ~ J " J '^('"1' ■■■iUp\z)(f){ui, ...,Up\z)du\ ■ ■ ■ dup. (1.7) 

The complex variables zi,...,zn of (1.7) are related with the disorder parameters of the OSBAE . The 
vector valued function (f){ui, ...,Up\z) is the quasi-classical limit of the wave vector $(ui, ...,Up\z). The 
Bethe ansatz for the Gaudin model was derived for any simple Lie algebra by Reshetikhin and Varchenko 
[20]. The scalar function X{ui, Up\z) determines the monodromy of ^'(^;i, . . . , Zn) and it is constructed 
from the quasi- classical limit of the A('u = Zk; u\, Up\z) and Ta{ui, - ■■ , Up\z) functions. This represen- 
tation of the A''-point correlation function shows a deep connection between the inhomogeneous vertex 
models and the WZNW theory. 
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In [21] this idea has boon apphcd on the periodic XYZ Gaudin model [22] in order to derive and solve 
the Knizhnik-Zamolodchikov-Bernard equation for S'?7(2)-WZWN model on the torus. [23]. In a series of 

(2) 

papers [24, 25, 26] we have considered the 19-vertices models, the osp{l\2) model and the twisted A2 
or Izergin-Korepin model and the s?(2|l)(^) model, respectively. In these works, the OSBAE procedure 
was used in order to find and solve the corresponding Gaudin models, as well as, the corresponding KZ 
equations. The algebraic structure of the periodic trigonometric osp{l\2) Gaudin model has been studied 
with details in [27]. 

Since the work of Sklyanin [28], the algebraic Bethe ansatz has been applied to various integrable 
models with non- trivial boundary conditions, which are specified by ii'-matrices satisfying the reflection 
equation and its dual [29] . The quasi- classical expansion of the corresponding double-row transfer matrices 
produces generalized Gaudin Hamiltonians with boundaries. 

Recently Gaudin models with non- trivial integrable boundaries have attracted much interest. Ini- 
tially the attention has been concentrated on Gaudin models with diagonal ii'-matrices [30, 31, 32, 33]. 
Nevertheless, in [34] the XXZ Gaudin model was solved with non-diagonal if -matrices and in [35] this 
result was generalized in order to solve the An-i Gaudin magnets with non-diagonal if -matrices. In this 
context, in [36] the XXZ model is also used to solve the Gaudin magnets with non-diagonal impurities. 

In this paper we present a systematic study of the algebraic structure of the 19-vcrtcx model based 
on the orthosymplctic Lie algebra osp{l\2), through the algebraic Bethe ansatz with diagonal if-matrices 
and we also show how the quasi-classical limit procedure gives us the corresponding generalized Gaudin 
magnets. Moreover, we will construct and solve the trigonometric osp(l|2) KZ equations. 

The paper is organized as follows. In Section 2, we present the osp(l|2) vertex model and its three 
diagonal if-matrices. In Section 3, the inhomogeneous algebraic Bethe ansatz with boundary terms is 
presented with details in order to derive the boundary off-shell Bethe ansatz equation for this vertex 
model. In Section 4 , the structure of the generalized Gaudin Hamiltonians is presented. In Section 5, 
we describe the corresponding off-shell Gaudin equations taking into account the quasi- classical limit of 
the results of the Section 3. In Section 6, data of the boundary off-shell Gaudin equation are used to 
construct solutions of the trigonometric KZ equation for the three cases. Conclusions are reserved for 
Section 7. Finally, commutation relations and some definitions are presented in appendices. 

2 The model 

Let V = Vo®Vi be a Z2-graded vector space where and 1 denote the even and odd parts respectively. 

s 

The multiplication rules in the graded tensor product space V differ from the ordinary ones by the 

s s 

appearance of additional signs. The components of a linear operator B e V result in matrix 
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elements of the form 

{A ® ByJ^ = (-I)f(/3)(P(«)+P(7)) A„^B0s. (2.1) 

S S 

The action of the graded permutation operator V on the vector |a) (H) |/?) G y y is defined by 

V \a) ® 1/3) = (-)f(")f(/5) 1/3) 4 |a) ^ (P);^'^ = (-l)f(")f(^)j„5 60^. (2.2) 
The graded transposition st and the graded trace str are defined by 

= (-l)^^^"' + '^^<''^^/3a, StrA = ^(-1)''(")A„„. (2.3) 

a 

where p{a) = 1 (0) if \a) is an odd (even) element. 
For the graded case the Yang-Baxter equation 

7^l2 (u)7^l3 {u + v)n23 (v) = n^z {v)ni3 {u + v)ni2 («) (2.4) 

and the reflection equation [28, 37] 

7^12(M - v)K:[{u)n2i{u + v)K^{v) = K^{v)ni2{u + v)K:[ {u)'R2i{u - v) (2.5) 

remain the same as in the non-graded cases and we only need to change the usual tensor product to the 
graded tensor product. 

In general, the dual reflection equation which depends on the unitarity and cross-unitarity relations 
of the 7^-matrix takes different forms for different models. For the model considered in this paper, we 
write the graded dual reflection equation in the form introduced in [38] : 

n'2{''*' i-u + v){K+y''^ (w)Mf l7^f2^"*^ {-u-v- 2p)Mi{K+y^^ (v) 

= {K+y''{v)MiTZll''*^-u -V- 2p)Mf \i^+)^*l(w)7^^*l^*=^(-u + v), (2.6) 

Now, using the relations 

7^12' (w) =-^1^21 (u)/i, 7^fl^''*'(M) ==/ii?i2(w)/i and IK+{u)I ^ K+{u) (2.7) 

with I = diag(l, —1, 1) and the property [M1M2, 7?.(w)] = we can see that the usual isomorphism [39] 

K-{u) K+{u) = K-{-u - py*M. (2.8) 

holds with the BFB grading. Here st^ denotes super-transposition in the space i. 

A quantum-integrable system is characterized by the monodromy matrix T(u) satisfying the funda- 
mental relation 

R{u - v) [T{u) T{v)] = [T{v) ® T{u)] R{u - v) (2.9) 
4 



where the intertwining _R- matrix is given by R{u) — VTZ{u). 

In the framework of the quantum inverse scattering method [15], the simplest monodromies have 
become known as C operators, the Lax operators, here defined by Caq{u) = TZaq{u), where the sub- 
script a represents the auxiUary (horizontal) space, and q represents the quantum (vertical) space. The 
monodromy matrix T{u) is defined as the matrix product of N Lax operators on all sites of the lattice. 



T(u) = jCaN{u)jCaN-l{'u) ■ ■ ■ j0.al{u). 



(2.10) 



The main result for open boundaries integrability is: if the boundary equations are satisfied, then the 
Sklyanin's transfer matrix [28] 



t{u) = str„ {K+{u)T{u)K-{u)T-\-u)) 
forms a commuting collection of operators in the quantum space 

[t{u),t{v)]=0, yu,v 



(2.11) 



(2.12) 



The commutativity of t{u) can be proved by using the unitarity and crossing-unitarity relations, the 
reflection equation and the dual reflection equation. In particular, it implies the integrability of an open 
quantum spin chain whose Hamiltonian (with -fl'~(0) = 1) is given by [28] 

k=i ^ "-^ «=o 

where the two-site terms are given by 



stroj^o+(0)gjv,o 
sivK+{Q) ' 



(2.13) 



Hk,k+1 = -l-Rk,k+l{u) 

du 



(2.14) 



•u=0 



in the standard fashion. 

The trigonometric solution of the graded Yang-Baxter equation (2.4) corresponding to osp(l|2) in 
the fundamental representation has the form 



n{u) 



1 

X2 



/ 
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Xl ) 



(2.15) 
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with non-zero entries [40]; 

Xi{u) = sinh(u + 2r]) sinh(u + 877) , X2 (w) = sinh u sinh(u + 877) 
Xsiu) = sinh u sinh('U + r?) , X4,{u) = sinh u sinh(u + 877) — sinh 2?7 sinh 3r] 
X5 (u) = e~" sinh 2r] sinh(u + 877) , 2/5 (u) = e" sinh 2r] sinh(u + 877) 
xq{u) = —ee~"~^'' sinh 2?7 sinh u, t/g (^t) = ee""'"^'' sinh 2?? sinh u 
xt{u) = e^" sinh 2?7 (sinh(M + 377) + c^'' sinh u) 

y^{u) = e" sinh 2?7(sinh(u + 377) +e'' sinh u) (2.16) 
satisfying the properties 

regularity : Ui^iO) = 1(0)'^^ P12, 
unitarity : ni2{u)nll''*^ {-u) = f{u), 
PT - symmetry : Pi2'R-i2{u)Pi2 = Kl''*' {u), 
crossing - symmetry : 7^l2(u) = UiTll*^ {-u - p)U^'^, (2.17) 

where e = ±1, f{u)= i^^(^) ^2(-m) ^^"^ ^ crossing parameter and U determines the crossing matrix 

M = U*U = M*. Here we have assumed that the grading of threefold space is p(l) = p{3) = and 
p(2) = 1, the BFB grading and we will choose the solution with e = 1. Moreover, we have verified 
that the three different sets of equations associated with the gradings BFB, FBB and BBF are mutually 
equivalent. 

This 7^-matrix is regular and unitary. It is PT-symmetric and crossing-symmetric, with p = 3r] and 




M = 1 . (2.18) 



Diagonal solutions K~{u) for the refiection equation (2.5) have been obtained in [41]. It turns out 
that there are three solutions without free parameters, being K~{u) = 1, K~{u) = F+ and K~{u) = F~, 
with 

=pe-2«/(±)(u) 

1 I , (2.19) 



Te^"/(±)(u) 



where we have defined 



By the automorphism (2.8), three solutions K^{u) follow as K~^{u) = M, K~^{u) = G"*" and K^{u) = G~ , 
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with 



=pe2«+%(±)(u) 
= { 1 I , (2.21) 

Te-2"-^''5(=^H") 



where we have defined 



„(+) = ^^Hu+lv) (_) _ cosh(u+|r?) 

We have thus nine possibihties for the commuting transfer matrix (2.11). We will only consider 
three types of boundary solutions, one for each pair {K~{u), K~^{u)) defined by the automorphism (2.8): 
(1,M), (i^+,G+)and {F-,G-). 

3 The boundary off-shell Bethe equation 

The boundary algebraic Bethc ansatz for a W-sitc inhomogcneoiis do\ible-row transfer matrix is obtained 
from the homogeneous ones by a local shifting of the spectral parameter u in the monodromy matrix 
T{u) and in its reflection T^^{—u): 

3 ^ s s 

Tia{u\z)= ^kii"^- ^N,V)^^k!k2^\u- ZN-l,V)^---^^kLiai'^- ^-^^'n) 

fei,...,feN-i = l (3.1) 

and 

3 

T,-\-u\z) = Yl + ^) ® 41 + ^2, ^) ^ • • • ^ 4^i, (« + ZN, V)- 

fci,...,fejv_i=i (3.2) 

where 

<f (« + ^„,r?) = ] rW(-^-^„,-^), (3.3) 

are the 3x3 matrices (obtained from the identification Caq{u) = TZaq{u) ) acting on the nth site of the 
lattice. 

The double-row monodromy matrix U {u) can be written as a 3 x 3 matrix 



/ Uii{u\z) Ui2{u\z) Ui3{u\z) 

U{u\z)=T{u\z)K-{u)T-\-u\z)= \ U2i(u\z) U22{u\z) U^sMz) \ (3.4) 

\ U^iiulz) U32{U\Z) U3:i{u\z) 

and it is more convenient to introduce new operators: 
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'Di{u\z) = Un{u\z), 

'D2{u\z) = U22{u\z) - fi{u)'Di{u\z), 

'D3{u\z) = U33{u\z) - f2{u)'Di{u\z) - f3{u)'D2{u\z), 

Bi{u\z) = Uij{u\z) {i < j), Ci{u\z) = Uij{u\z) {i > j), i,j = 1,2,3 (3.5) 
The usual reference state is 

|0)=|^oj|)^oj®---|)^oj (3.6) 
which is a highest vector of W(u) 

/ Xi{u\z) |0) * * 

U{u)\0) = * I (3.7) 

V X3{u\z)\0) 

where the symbol (*) denotes states different from and |0) and 

N 



= fcn(n)n 



Xlju - Z„) Xl{u + Zg) 1 
^J-^ X2iu - Za) X2{U + Zg) f {u + Zg) '' 



^ 1 

X2{u\z) = [k22{u) - A;ii(u)/i(u)] TT 

a=l + 



)' 



Mu\^) = [k,3iu)-k,2iu)f3iu)-K,iu)Mum^^ (3.8) 



where 



r f ^ ^ 1/5 (2m) f m = ^1(2^)2/5(2^) - X5{2u)y7{2u) 

~ xi{2u)' ■^3W- xi{2u)x4i2u) + X5i2u)y5i2u)' 

f (,A = ^7 (2m) , , ^ Xii2u)y'r{2u) + y1{2u) 

' xi{2uy ■'^^ ' xi{2u)xi{2u) + X5{2u)y5i2u)' ^ ' 

Using the K~^{u) reflection matrices (2.21) and the BFB grading, the double-row transfer matrix 
(2.11) has the form 

t{u\z) = k^^{u)Un{u\z) - k}2{'^W22{u\z) + k^^{u)U33{u\z) 

= ni{u)Vi{u\z) + n2{u)'D2{u\z) + Q.3{u)V3{u\z) (3.10) 

where 

fli{u) = ktiiu) - fi{u)k^2iu) + f2{u)k^3{u), 
^2{u) = -k^2iu) + f3iu)k^3{u), 

^s{u) = k+{u). (3.11) 



8 



The monodromy matrix lA [u) satisfies the relation 

7^l2(u - v)Ui{u)'Ti2i{u + v)U2{v) = U2{v)ni2{u + v)Ui{u)Tl2i{u - v) (3.12) 

from which wc obtain the commutation relations for its entries. The commutation relations which par- 
ticipate effectively in the Bethe ansatz construction are presented in the appendix A. 

The n-particle state is defined as an operator valued function through the recurrence relation 



^n{u,... ,Un\z) = Bi{ui\z)<bn-l{U2,--- ,Un\z) 

n 

+B2{ui\z)^Xl{Ui\z)Fi\ui,... ,M„)*„_2(W2, • • • ,Ui,... ,Un\z) 



i=2 

-W 

i=2 



+B2{ui\z) ^ X2(Ui\z)F2\ui, ■ ■ ■ ,Un)'^n-2{U2, ■ ■ ■ \z) (3.13) 



with = |0) and ^'i(ui) = B{ui\z) |0). Here Ui denotes that the rapidity ttj is absent. 

It was shown in [42] that this operator is normal ordered satisfying n — 1 exchange conditions 

*„(Ul, ... ,Ui, Uj+i, . . . , Un\z) = U){Ui, Wi+i)4'„(ui, . . . , Uj+i, U,, . . . , Un\z) (3.14) 

provided the functions Fa\ui, . . . ,Un) are given by 

i—l n 

Fli\ui,... ,Un) = '[\.Oj{Uj,Ui) Yi aal{Ui,Uk)Gd„{ui,Ui), (« = 1, 2) (3.15) 

j=2 k=2,kjti 

where 

X3{u - v)x4,{u - v) - xe{u - v)ye{u - v) 
w{u, v) = ^ r — --^ 

Xl{U — V)X3{U — V) 

ui{u2,v)ui{u,v) = 1 (3.16) 

and 

xe{u - v) X2i2v) xe{u + v) 

X3{U — V) Xi[2v) X2{U + V) 

Using the commutation relations listed in the appendix A, one can find the action of the operators 
Va{u\z) , a = 1,2,3 , on the n-particle state. In order to do so we can start with the one-particle state, 
then the two-particle state and so on. This procedure is very extensive and tedious but the final result 
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can be presented in a compact form 

n 

T>a{u\z)'^n{ui, . . . ,Un\z) = Xa{u\z)Y[aal{u,Ui)^n{ui, ■ ■ ■ ,Un\z) 



i=l 

n i—1 I n 

+ ^ Xi{ui\z)aa2{u,Ui)Y\_an{ui,Uk) 

i=l j=l \ kjti 

\ 

+ X2{Ui\z)aa3{u,Ui)Y\_a2l{Ui,Uk) 1 Bi{u\z)'l!n-i{ui\z) 

k+i ) 

n i—1 I n 

+ (1 - ^ Yl. ^("i' Xi{Ut\z)aa4{u, Ui) aii{Ui, Uk) 

+ X2{Ui\z)aa5{u,Ui)Y\_a2l(Ui,Uk) I B2,{u\z)^ n-l{ui\z) 
k^i J 

n—ln\ n n 

+ ^ ^ < Xi{ui\z)Xi{uj\z) Yl aii{ui,Uk) Y\ aii{uj,Ui)Hai{ui,Uj) 

i=l \^ k^ij 

n n 

+X2{ui\z)Xi{uj\z) Y[ a2i{ui,Uk) Y\_ an{uj,ui)Ha2{ui,Uj) 

n n 

+X-i{ui\z)X2{uj\z) aii{ui,Uk) a2i{uj,ui)Has{ui,Uj) 

n n 
+ X2{Ui\z)X2{Uj\z) a2l{Ui,Uk) W a2l{Uj,Ul)Hai{Ui,Uj) 

xJJa;(Mfc,Uj) JJ co{ui,Uj)B2{u\z)'bn-2{ui,Uj\z) (3.18) 

fe=l l=l=ii 

The action of the transfer matrix t{u\z) on this n-particle state defines the so-called off-shell Bethe 
ansatz equation [18]. In our case it has the following form 

n n J — 1 

t{u\z)^n{Ul, . . . ,Un\z) = AMUU . . . ,n„\z) + E-^i""'^*(n-l) + EE^?"'^<-2) 

3 = 1 j=2 1=1 (3.19) 

Let us explain a bit more the right hand side terms of this equation: In the first term the Bethe vectors 
(3.13) are multiplied by c-numbers A„ = A„(u, ui, ...,Un\z) given by 

3 n 
An = ^Q.a{u)Xa{u\z) JJ aai(M, Ufe). (3.20) 

a=l k=l 

The second term is a sum of new operator valued functions 

j-i / 3 3 \ 

*(n-l) = H'^^^'^'^j) ^^a{u)aa2{u,Uj)Bi{u\z) + ^na{u)aa4{u,Uj)B3{u\z) *(„_i)(«j) 

fe=l \a=l a=2 / (3.21) 
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multiplied by scalar functions 

n n 

^ Xi{uj\z) JJ an{uj,Uk) + Q{uj)X2{uj\z) a2i{uj,Uk) (3.22) 

k^j k^j 

where we have defined 



^a{u)aa3{u,Uj) ^ na{u)aa5{u, Uj) 

&{uj) = ^ = ^ . (3.23) 



Finally, the last term is a coupled sum of a third type of operator valued functions 

*(l-2) = n '^("fc''"') n '^("fc''"j)^2(w|^)*(n-2) {ui,Uj) (3.24) 
k=l k=l,^l 

with intricate coeSicients 

n n 3 

jr(n-2) ^ X,(ui\z)Muj\z)l[ 

kjil k^,j a=l 

n 71 3 

+X2{ui\z)Xi{Uj\z)Y\_a2l{ui,Uk) JJ aix{Uj,Uk) ^ 9.a{u)Ha2{ui,Uj) 
kjtl k^j a=l 

n n S 

+Xi{ui\z)X2{uj\z)Y[an{ui,Uk) JJ a2i{uj,Uk) ^ Vla{u)Has{uuUj) 

kjtl k^j o;=l 

n n S 

+X2{ui\z)X2{uj\z)Yla2i{ui,Uk) Yl a2i{uj,Uk) ^ na{u)Ha4{ui,Uj) (3.25) 

kjtl kjtj a=l 

The functions aij{u,v) are amplitudes which came from the commutation relations (appendix A), 
while Hoij{up, Uq) {a = 1, 2, 3, j = 1, 2, 3, 4) are cumbersome functions of this amplitudes which we have 
left to the appendix B. 

In the usual algebraic Bethe ansatz method, the next step consist in impose the vanishing of these 
unwanted terms in order to get the eigenvalue problem for the inhomogeneous transfer matrix t{u\z): 
We impose = 0, which implies = 0. From these elimination conditions wo can see that 

\l/„(tti,... ,Un\z) is an eigenstate of t{u\z) with eigenvalue (3.20) provided that the rapidities Uk are 
solutions of the inhomogeneous Bethe equations 

$^=-e(..) n ^^4^, ik^l,2,...,n) (3.26) 
where the factors @{uk) are easily identified with the phase shifts at the boundaries. 
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4 The Gaudin Hamiltonians 



We recall that the osp{l\2) has three even (bosonic) generators H, and two odd (fermionic) generators 
V"^ , whose non- vanishing commutation relations in the Cartan-Weyl basis reads as 

= ±X±, [X+,X-] = 2H, 

{y±,i/±} = ±ix±, {v+,v-} = -Ih. (4.1) 

The quadratic Casimir operator is 

C2 = H^ + l{X+,X-} + [V+,V-], (4.2) 

where {• , •} denotes the anticommutator and [• , •] the commutator. 

The irreducible finite-dimensional representations pj with the highest weight vector are parametrized 
by half-integer s = j/2 or by the integer j e N. Their dimension is dim{pj) = 2j+l and the corresponding 
eigenvalue of the quadratic Casimir is j{j + l)/4 = s{s + 1/2). 

The fundamental representation has s = 1/2 and its matrix realization in the BFB grading is 

H = 






f 
























/ 



= - 1 , y- = - -1 . (4.3) 



In this section we will consider the theory of the Gaudin model with boundary terms. To do this we 
need to calculate the quasi-classical limit of the results presented in the previous section. 

We can begin expanding (up to an appropriate order in rj), the diagonal Lax operator entries of the 
monodromy matrix T{u\z) 



= In + 2rj coth(« - zn)nn + 2rf (nl + I ^" , ) + 0(7,3) 



6^ = Ir.-2rf ^^",, +o(r?3) 
smh(M — Zn) 



4^) = T„-2r? ootHu-z^)Hn + 2rf{K + \ ^ '^\, ) + o{ri% (4.4) 
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For the elements out of the diagonal we have considered a different order in ry 



,(n) 
-12 



-23 



-2tj 



2V ^ 



sinh(u — Zn) 

„-«+z„ 



21 



277 



sinh(u — z„) 



V+ + o(r?2), 



(n) 
13 



27? 



sinh(w — 



V-+o(r?2), 



-32 



2»? - 



sinh(M — Zn) 



V++o(r?2), 



(") 

31 



2r? 



sinh(u — Zn) 



(4.5) 



sinh(u — Zn) 

where V± = 2F±, Af^ = 2X± and W = 2H. 

The corresponding expressions for the reflected monodromy matrix T~'^[—u\z) can be read from (4.4) 
and (4.5) by the substitution —r\ and u — > —u. For sake of simplicity we can remove the global factor 
by considering its expansion as 11(1//) w 1 + o(?7^) instead of a normalization. 

The quasi-classical expansions for the entries of the K matrices (2.18)-(2.22) will be written in the 
form 



fcaaW = W + ^^ai'^M + ^^'^^1'^^) + o(??'), (« = 1, 2, 3). 



(4.6) 



Substituting these expansions into the matrix operator V = K^{u)U(u\z) we will get a quasi-classical 
expansion for its diagonal entries: 



VMz) = 1 + 1riV^2{u\z) + 4r?2yj2)(^|^) + 0(^3-, 
where the first order terms are very simple 

N N 
+ ^2 COth(u - Za)'Ha + ^ COth(u -|- Zi,)'Hb, 
0=1 6=1 

N N 
— ^ COth('U - Za)'Ha — ^ COth('U + Zb)Hb 



(4.7) 



Vll\u\z) 



Vs^P{u\z) 



(4.8) 



a=l 



6=1 
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and the second order terms are very complicated 

+ ^ |coth(u — Za) coth(u — Zb)Ha ^Hb + coth(u + Zb) coth(u + Za)'Hb (8) I 



+ \ktl'\u)k;l'\u) 



a<b 



a<5 



sinh(u — Za) sinh(w — Zb) 

,2 N N 



sinh(u + Zb) sinh(u + z^) 



1-7/2,3 Til —Hb 



4 sinh(M + ZbY 



1 

+ 2 



6=1 



kti°\u)k^i''\u) + ktl^\u)k^l:''\u)\ <j ^ coth(« - + ^ coth(u + Zb)nb 

N N 

+ ^ coth(u - 2;o)'Ha ^ coth(u + 2;fe)W6 



a=l 



b=l 



a=l ^ ' 6=1 ^ ' 



N 



AT 



^ sinh(u - Za) ^ sinh(u + Zb) 



(4.9) 



1 [fc2+^^)(t.)fc22^°H«) +4'°'(«)fc2"2^')(^.) 



AT 



a<5 



sinh(w — Za) sinh(u — Zb) 2 sinh(w — Za)"^ 



JV 



a<6 



sinh(u + Zb) sinh(w + Zb) 2 sinh(u + Zb)^ 



^ sinh(u - Za) ^ sinh(u + ^b) 



JV 



V- 



N 



^_U-Zb^ + 



^ sinh(u - Za) ^ sinh(u + Zb) ' 



(4.10) 
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■ ^ |coth(u — Za) coth(u — Zb)Ha ^Hb + coth(u + Zb) coth(u + Za)'Hb (8) I 



a<b 



a<b 
N 



a=l 
1 



sinh(M — Za) sinh(M — Zb) 
,2 . ^. N N 



sinh(u + Zb) sinh(u + Zg) 



1^/2 3 Ha+Ha 



2 " 4sinh(w-^a)V ' ^ " 4 sinh(M + 2^)2 

( N N " 

ktP{u)k^P{y') + ktP{u)k^P{u)\ I J2 - Za)na + ^ COth(u + ^b)^^ 



AT 



+ ^ COth(u - Za)Ha ^ COth(u + Zb)Hb 



a=l 



b=l 

Af I Af I 



+*3r'(»fe-r(»)EiiT7;7^E 



^ sinh(w - Za) f—f sinh(w + z^) 

a=l 0—1 ^ 



+kt}°\u)k^l'\u) V ^ E - 



^ sinh(u - Za) ^ sinh(u + Zb) 



(4.11) 



Here we have used the identities kaa\u)kaa\u) = 1, a = 1,2,3., which hold for the three boundary 
solution due to the isomorphism (2.8). 

Now, the quasi-classical expansion for the double-row transfer matrix (2.11) has the form 



t{u\z) = Vll{u\z) - V22{u\z) + V33{u\z) 

= 1 +87?V^)(w|z) + o(r/3). 



(4.12) 



where we notice that the first order terms in r] are canceled. 

The residue procedure of t^'^\u\z) at the point u = Za results in the following generalized Gaudin 
Hamiltonians 
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Ga — - (j^ii \^a)kii \za) + \za)kii \za) k^^ \za)k^^ \za) k^^ ^(-20)^33^ \za)^ Tia 

^ 1 ( 1 

^ sinli(2;a - Zb) [_ z 
+ e"'-'-V+ I) Vb - e^^-^'-V" 1) Vfc+I 

1 ^ 1 

+ E . w ^ , I (fc^"r(^a)fcrr(^a) + fear's °^(^a)fc3"3^°^ (-2a)) COSh(^„ + ^ Wfe 

2 f-^ SmIliZa + Zb) 1- V / 

0—1 



- + k;l'\z,)k-,^'\z,)) e-^-^'V- ® V+ 

+ (^k+,^°\z,)k-l°\z,) + k+,^°\za)k-,^'\za)) e^»+^^V+ ^ V,-} (4.13) 

which satisfy the following integrable conditions 

EG„ = 0, ^ = ^, [G„,G5]=0, Va,6. (4.14) 

a— 1 

It means that wo have derived N — 1 independents non-local integrable Hamiltonians for each of the 
three boundary solutions (1,M), (F",G") and {F+,G+). 

Let us look at the right hand side term of (4.13): The first term represents the diagonal boundary 
term, the first sum is nothing else but the sum of classical r-matrices Ylb=ia ^ab{zz — Zb), as it could be 
expected and the second sum is also a sum of classical r matrices but with the boundary effects. This 
difference emerges due to the breaking the symmetry of the double-row structure when we choose the 
point of a particular row to take the residue. Therefore, we have the following generalized form 

Ga = (b.t.)„ + rab{Zz -Zb) + Yl ^'abiZz + Zb) (4.15) 

b^a b^a 

where (b.t.)a is the boundary term at the site a (with residue at u = Za), rab{zz — Zb) is the "bare" 
classical r-matrix and '"'jj/j^z + Zb) is the classical r-matrices "dressed " with the boundaries. 

In the next section we will use the data of the boundary algebraic Bethe ansatz in order to find the 
exact spectrum and eigenvectors for these Hamiltonians. 
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Next, let us write the entries of the monodromy matrix U{u\z), up to their classical contribution: 



N 
a=l 






\ sinh(u 


- Za) 


N / 

= 2r?E 

a=l \ 


fc33^"''(M)e" 


-U+Za 


sinh(u — 


Za) 


N 1 

= 2.E 

a=l \ 




-U-\-Za 


sinh(u — 


Za) 



(4.16) 



for the creation operators and 

C^{u\z) = 2r?E f + fe22^%)e-"-^° A V+ +o(r?^) = 2r,V+ {u\z) + o{r,') 

^V ^^^^"~'^"^ smh(M + 2:„) J " ^'^ /VI/ \i J 



N 



sinh(u — Za) 


sinh(M + Za) 






sinh(u — Za) 


sinh(« + Za) 


k^l°\u)e--''^ 


kss^'Hu)e—^^ 



(4.17) 

for the annihilation operators. For the diagonal entries of U{u\z) we have: 



V3{u\z) = k^P{u) + 2'nn{-u\z) + o{rf) (4.18) 



where 



'H{u\z) = k^l°\u) E (cosh(u - Za) + cosh(M + Za)na + ^A;-/^^(u). (4.19) 

0=1 

Here we have used the relations A;^^"^ (u) = fc^/°^(— u), fc^^°^(u) = k^2^\—u) and fc^^^^(w) = — fc^/^^(— u). 

Note that with these expression we are defining the generators for a "double" osp{l\2) Lie superal- 
gcbra. The defining relations of this algebra can be obtained from the quasi-classical limit of (3.12) or, 
equivalently, from the quasi-classical limit of the commutation relations of the appendix A. 

5 The off-shell Gaudin equation 

The off-shell Gaudin equation is defined from the residue of the quasi-classical limit of the off-shell Bethe 
ansatz equation (3.19). In this limit we first consider the quasi-classical expansions of the states (3.13), 
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(3.21) and (3.24) in terms of the operators (4.16) 

*„(ui, Un\z) = (-2r?)"$„(ui, Un\z) + o«+i), (5.1) 



= 2(-2,7)"+H-l)^'|^[fciT^(«) + fc+(°)(n)fc+^°^(«)]- 



sinh(u — Uj) 

+ l[42^°\u) + °^ V-{u\z)^n-i{uj\z) + o(r,"+2) (5.2) 

= {-2vr-H--^y^''^-m^n-2{ui,uj\z) + o{vn (5.3) 
where the quasi-classical Bethe state $n(ui, ...,u„|2) is given by the recurrence formula 

^n{ui,...,Un\z) = V~{ui\z)^n-l{U2,--,Un\z) 

-^-(-^i-)E 3U(.,-.,) ^" ^"-^^"-^^"-'^^ 

''^ / \jp — Ui—Uj 

E sinh(., +.,) '- ^"--^^"-^^"-'^^ ^'-'^ 

with $0 = |0) and ^i{ui\z) = V~{ui\z)^o- 

Next we consider the quasi-classical expansions of the c-numbers functions (3.20), (3.22) and (3.25) 

A„ = l + 2(-2r?)2Al2)+o(r;3) 

= 2{-2vr hi+o{n% (5.5) 
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where 

n 1 



+ ktl'Hu)k-l°\u) - kU'\u)k-l°\u) + kti'\u)k-,^'\u)} 



1 (fc+/")(z.)fcr/''(«) + kt}'\u)k;l"\u) - k+^'\u)k-,^'\u) - kt,^'\u)k-,^'\u) 



+ 

4 

' AT 



X 



^^[C0th(u — Za) + COth(« + Za)] — ^^[coth(u — Uj) + COth(u + Uj)] 
ya=l j=l 



V-/w ^3 1 3 1 
+ } [ COth(u - ZajCOth{u + Za) + -^-r-, T? + 7^T7 i 

^ V 4 sinh(u - Za)^ 4 sinh(w + Za)^ 

" / 1 1 , 

+ ( coth(u — Uj) coth(u + Uj) + - coth(u — Uj)"^ + - coth(u + Uj)' 



3=1 

N 



+ {C0th.{u — Za) COth.{u + Za)){coth.{u — Zb) COth.{u + Zb)) 

a<5=l 
n 

+ {coth{u — Ui) coth{u + Ui)){coth{u — Uj) coth{u + Uj)) 
i <j=i 

N n 
— ^^(C0th(u — Za) + COth(u + Za)) ^^(coth(u — Uj) + COth(u + Uj)) (5.6) 
a=l j=l 

N n 
fj = - ^(coth(Uj - Za) + COth(Uj + Za)) + ^(coth(Uj - Uk) + COth(Uj + Uk)) 

a=l k^j 

Note that in these relations we have expanded (3.23) as 9{u) = 6'(''^(u) + rjO^-^^u) + oijf') and used the 
identity k-^i^\uj) + k22^\uj)9^^\uj) = which holds for the three cases. 

The contribution of involves the expansions of the complicate functions Haj presented in the 

appendix B. After some algebraic manipulations we have 

^ sinh(2M) [ sinh(2Mj) fcn°VO + fe^2°^(Mi) n,-(o)/ wi 

sinh(uj + ui) sinh(ui - Uj) | sinh(u - ui) sinh(M + ui) fc+(o) (y^.) + fc+(o) (y^.) ' ' 

^^^^^i^^^^^ fen"H'0) + 4^'''(".) [fe-W(^^.)^^.]l (5_8) 



sinh(u - «j) sinh(u + u^) fc^"/°^(u() + k'^^^\ui) 



19 



Substituting these expressions into the (3.19) and comparing the eoefhcients of the terms 2(— 2ry)""'"^ 
we get the first non-trivial consequence for the quasi-classical limit of the OSBAE: 

j^i smn^u Uj) 



Here we have used the identity 

sinh(2t;) 



(5.10) 

Note that in (5.9) the contributions from and *('^_2) were combined in order to define two new 



sinh(u — v) sinh(u + v) sinh(u — v) sinh(u + v) 
'(n-i) and 



vector valued functions 

ei = i[fc+r(u)+^-r^(«)fc+^"^(n)]V-(z^|z) $„_!(«,>) 



(5.11) 



and 



©^2 = ^[fc2'2^°)(t.) + fciTH«)4^°H«)]V-(t.|^)*„-l(«,|^) 



' ' smhMfc-Mj smhK + Wj A;+/°^(ufe) + A;+^°^(wfe) 



fe=i 



(5.12) 



where fc = A: -I- 1 for k < j and k — k for k > j. 

Finally, we have the off-shell Gaudin equation taking the residue of (3.19) at the point u — Za- 



Ga^n{ui,...,Un\z) = ga^n{ui, ...,U„\z) + ^{-1)' t[ k-^}°\uj) f j](f>[ 



+ E(-l)^-K(^[^n^"KO/.^ 
1,2, ...,iV (5.13) 
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(2) 

where ga is the residue of A„ at the point u = Za 



+ 4 



sinh(22;a) 

(^1/ ^(•^0)^1/ ^(•^o) + ^11 ^(-20)^1/ ■'(•^^o) ~ ^33^ ^(-20)^33'' \za) — \za)k33 \za)^ | 

N n 
+ COth(22;a) + ^(C0th(2a - Zb) C0th(2:a + Zb)) - ^(C0th(2:a - Uj) + C0th(2o + Uj)), (5.14) 

The new functions </'^2) are given by 4>i^2) ~ ^^^u=za^i(^2)' 

In this way we are arriving to the main resuh of this paper. The equation (5.13) governs the Gaudin 
theory: From its reduction to an eigenvalue problem we can find the eigenvalues and corresponding 
eigenvectors for A'' — 1 commuting Hamiltonians Ga (4.13). It means that Qa is the eigenvalue of Ga 
with eigenfunction $„(ui, u„) provided ui are solutions of the equations fj = . Moreover, as we will 
see in the next section, the off-shell Gaudin equation defines the conditions on the monodromy of the 
correlations functions in order to be solutions of the differential equations known as KZ equations. 

6 The Knizhnik-Zamolodchickov equation 

The Knizhnik-Zamolodchikov (KZ) differential equation 

^9|(£) ^ (g 1) 

dZi 

appeared first as a holonomic system of differential equations on conformal blocks in a WZW model 
of conformal field theory. Here ^{z) is a function with values in the tensor product Fi • • • (g) V)v of 
representations of a simple Lie algebra, k = k + g , where k, is the central charge of the model, and g is 
the dual Coxeter number of the simple Lie algebra. 

One of the remarkable properties of the KZ system is that the coefficient functions Gi{z) commute 
and that the form u = J2i Gi{z)dzi is closed [20]: 

Indeed, it was indicated in [20] that the equations (4.2) are not just a flatness condition for the form ui 
but that the KZ connection is actually a commutative family of connections. 

In this section we will identify Gi with the bounded osp(l|2) Gaudin Hamiltonians Ga presented in 
the previous section and show that the corresponding differential equations (6.1) can be solved via the 
boundary off-shell Bethe ansatz method. 

Let us now define the vector- valued function \E'(zi, zat) through multiple contour integrals of the 
vectors (5.4) 

^{zi, Zn) = ■ " ^{u\z)^n{u\z)dui...dun, (6.3) 
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where X {u\z) = X (ui, tin, z\, zm) is a scalar function which in this stage is stih undefined. 
We assume that \E'(^i, z^v) is a solution of the equations 

«^^%^^ = Ga*(^i,...,^iv), a =1,2,. ..,7V (6.4) 

where Ga are the Gaudin Hamiltonians (4.13) and k is a constant. 
Substituting (6.3) into (6.4) we have 

where we are using a compact notation f {o} du = § . . .§ {o} du\ ■ ■ ■ dUn- 

Using the quasi-classical limit of the commutation relations one can derive the following non-trivial 
identity 

dza ^ au; Vsinh(ui - ^;„) smh.{ui + Za)) 
which allows us write (6.5) in the form 

dZa J y aza ^ smh(2:a-M;) smh(za + u/) J 

/ f '^(«l'^)[^^7— ^ + ^^T^^^]) '^"^ (6-7) 
^ J dui\ ^ ' ^^smh(^;a -wi) smh(^;a «;) V 

It is evident that the last term of (6.7) is vanishes, because the contours are closed. Moreover, if the 
scalar function A'(ulz) satisfies the following differential equations 

we are recovering the off-shell Gaudin equation (5.13) from the first term in (6.7). 

Next we can solve (6.8) in order to find the function X{u\z) which determines the monodromy of (6.3) 
as solution of the trigonometric KZ equation. 

Summarizing, let us consider the main results for the three boundary solutions (1, M), {F~^, G"*") and 
iF-,G-) 

6.1 The (1,M) Solution 

The data for this case are 

k;l'\u) = k-^'\u) = k-f \u) = 1, k;l'\u) ^ k-l'\u) = k-,^'\u) = 

(u) = (u) = fcg+f) {u) = l, 4(1) {u) = -ktl'^ {u) = 2, {u) = 

0"Hu) = -I, <,<■>(„) = -2j^ (6.9) 
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Substituting (6.9) into (4.13), (5.4), (5.14) and (5.7) we have the fohowing generahzed Gaudin Hamilto- 
nians 

^ 1 r 1 

+ e^^-^-X- ® X+) + e-^-+^»Vt ® n- - e^"-^^V- ® V+} 

^ 1 r 1 

+ E 3inh(.. + ..) I'^^^^^^" + ^ + 2 (^^""^^"" ^ 

a = 1,2,. ..,7V. (6.10) 
with eigenfunctions 

^n{u\,...,Un\z) = V~{Ui\z)^n-l{u2, ■■■,Un\z) 

-^"("^i^)Ei47^*-(%i^) 



(_)ie-«i-". A 

-^^ yui\z) 

and eigenvalues 



AT n 

5a = 2 coth(22:o) + ^ [coth(^;a - Zh) + coth(2:o + 2:5)] - ^[coth(^;(i - Uj) + coth(^;a + Uj)\ 

67^0=1 j=i (6.12) 

provided that 

TV n 

^^[C0th(uj — Zq) + COth(Mj ■ + ^a)] = [coth(Mj — Mfe) + COth(Mj + Mft)] (6.13) 
a=l /s#j=l 

The corresponding KZ solution is 

^(z\, ...,Zn) = f ■ ■ ■ f X{u\z)^n{u\z)dui...dun, (6-14) 



where 

N N 

X{u\z) = [sinh(2z„)]i/« [] Yl [sinh(z„-Z(,)sinh(z„ + 2b)]i/'= 

a=l 6=a+l 

r?, n N n 

X [sinh(wj - Uk) sinh(uj + Ufc)]^/" ]J[sinh(za - m^) sinh(2;a + Wj)]"^/" 

j=l/s=j+l 0=1 j=l 

(6.15) 
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is the solution of the equations (6.8). 

Here we note that the if-niatrix solution (1,M) is a quantum- algebra invariant solution. Comparing 
(6.10) with (4.15) we can see that r = r. It means that there is no boundary contributions on the bulk 
[43]. 

6.2 The (F±,G±) Solutions 

The data for the (i^+, G+) solution are 

fc-/^) (u) = -3e-2" coth u, k22^^^ (u) = 0, k^^^^'' {u) = -Se^"" coth u 

{u) = -46^" + Se^" coth u, fc^j^^^ (u) = 0, k^P (u) = 4e-2" + Se'^" coth u 
e^°^u) = e-2", ^(^)(u) =e-2"tanhu (6.16) 

while for the {F~ ,G~) solution we have. 

fcr/°^(n) = e-2«, k^l''\u) = 1, fcsa^^Hu) = e^" 

fc-/^) (u) = 3e-^" coth u, k^^^^ (w) = 0, k'^^^^ {u) = 3e^" coth u 

ktl'\u) = e^^ 42^°)(«) = 1, 4(°)(«)=e-2" 

(u) = 4e2" - 3e2" coth u, fc+*^^ (u) = 0, /t+g^^' (w) = -4e-2" - 3e-2" coth u 

0(o)(m) = e-2", 0(i)(m) =e-2"tanhM (6.17) 

We will present the summary of these two cases in a compact notation. 
Substituting (6.16) ((6.17)) into (4.13), (5.4), (5.14) and (5.7) we have 

^ 1 f 1 

Ga = 2Ha + ^ ^.^^^^^ ^ |cosh(^„ + z,)na k>n, + - (b-^^^ ^ 

+ e^^^e-^"-^" A*- I) + e (e-2^<'e^»+^''V+ I) V^- - e^^oe-^^-^^V" ® V^') } 
+ E 3inh(.'-..) {^°^^(^« - ^ + ^ (^"^""^^'^«' ^ 

+ e^o-^'-A'- ® + e-""+"''V+ ® Vfe- - e^-'^^V" ® V^+j (6.18) 
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with eigenfunctions 



j=2 

and eigenvalues 

N 



^ sinh(ui - Uj) 



9a = [C0th(2;a - Zj,) + COth.{Za + Zb)] - ^[coth(Za - Uj) + COth{Za 



+2 coth(22;a) + ^(1 + e) coth(2;a) + ^(1 - e) tanh(^;„) (6.20) 



provided that 



JV n 

^^[coth('Uj — + coth('Uj + ^a)] = [coth{uj — Uk) + coth{uj + Uk)] 

o=l kTtj=l 

_i(l + e)coth(u,) - i(l -e)tanh(«,) (6.21) 

The KZ solution is 

"beizi, zn) = j> ■■■ j> X{u\z)^n{'U'\z)du-i...dUn, (6.22) 

with 

X{u\z) = [cos\i{uj)]-^^-'^'^''[cos^Za)]^^-'^'^''[smh{uj)^^^^ 

N N n n 

X JJ JJ [sinh(^;a - Zb) sm]i{za + Zb)]^^'^ JJ JJ [sinh(uj - Uk) sm]i{uj + Uk)Y^'^ 

a=lb=a+l j=l fe=j+l 

N n 

X Yl JJ[sinh(z„ - Uj) smh{za + Uj)]-^/" (6.23) 

a=lj=l 

where e = 1 for {F~ ,G~) case and e = — 1 for {F~^, G+) case. 

7 Conclusion 

In this paper a detailed analysis of the boundary quantum inverse scattering method is applied in order to 
derive the off-shell Bethe ansatz equation (3.19) for a graded 19- vertex model based on the orthosymplectic 
Lie superalgebra osp{l\2). The boundaries ares given by three of solutions of the reflection equation and 
its dual, related in pairs by the isomorphism (2.8). The quasi-classical limit of this equation results in the 
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off-shell Gaudin equation (5.13) and wc have emphasized its importanee in the resolution of the Gaudin's 
theory. From this equation are following the energy spectrum of the Gaudin magnets and the monodromy 
function for the integral representation of the correlation function as solution of the KZ equation. 

There are several issues for which this paper could be useful: By the method presented in this paper 
one can include the other three 19- vertex models namely, the twisted s^(2|l)(^) model , the Zamolodchikov- 
Fateev model and the Izergin-Korepin model for which the algebraic Bethe ansatz with diagonal boundary 
are known [42, 45, 46]. We also would like to know how to connect the rational limit of our three 
trigonometric results presented in this work can be related with the osp{l\2) conformal field theory with 
conformal boundary conditions [44] and finally, should be really interesting to try extended these results 
in order to include all non-diagonal if-matrices classified in [41]. 

Acknowledgment: This work was supported in part by Fundagao dc Amparo a Pcsquisa do Estado 
de Sao Paulo-FAPESP-Brasil and by Conselho Nacional de Desenvolvimento-CNPq-Brasil. 

References 

[1] R.J. Baxter, Exactly solved models in statistical mechanics, Academic Press, London, 1982. 
[2] M. Gaudin, La fonction d'onde de Bethe, Masson, Paris, 1983. 

[3] V.E. Korepin , N.M. Bogoliubov, A.G. Izergin, Quantum Inverse Scattering Method and Correlation 
Functions, Cambridge- 1993. 

[4] E. Abdalla, M.C.B Abdalla, K. Rothc, Nonperturbative Methods in Two-dimensional Quantum Field 
Theory, 2nd Edition, World Scientific, Singapore, 2001. 

[5] M. A. Semenov-Tyan-Shanskii, Func. Anal. Appl. 17 (1983), 17. 

[6] A. A. Belavin, V.G. Drinfeld, Func. Anal. Appl. 16 (1982), 159. 

[7] M. Gaudin, J. de Physique 37 (1976), 1087. 

[8] M. C. Cambiaggio, A. M. F. Rivas and M. Saraceno, Nucl. Phys. A624 (1997), 157 [nucl-th/9708031] 

[9] L. Amico, Di Lorenzo and A. Osterloh, Phys. Rev. Lett. 86 (2001), 5759 [cond-mat/0105537]. 

[10] J. von Delft and F. Braun, Superconductivity in Ultrasmall Grains: Intriduction to Richardson's 
Exact Solution (1999) [cond-mat/9911058]. 

[11] J. Dukelsky, S. Pittel and G. Sierra, Rev. Mod. Phys. 76 (2004), 643 [nucl-th/0405011]. 



26 



[12] N. Seiberg and E. Witten, Nucl. Phys. B426 (1994), 19 [hep-th/9407087] . 

[13] V.G. Knizhnik, A.B. Zamolodchikov, Nucl. Phys. B247(1984), 83. 

[14] H.M. Babujian, R. Flume, Mod. Phys. Lett. A9 (1994), 2029 [hep-th/9310110]. 

[15] L.D. Faddeev, L.A. Takhtajan, Uspekhi Mat. Nauk 34 (1979), 13. 

[16] W. Galleas and M. J. Martins, Nucl. Phys. B699 (2004), 455 [nlin.SI/0406003]. 

[17] A. Lima-Santos, J.Stat.Mech. 0607 (2006), POOS [nlin.SI/0602003]. 

[18] H.M. Babujian, J. Phys. A: Math. Gen26 (1993), 6981 [hep-th/9307062]. 

[19] V. Schechtman, A. Varchenko, Invent. Math. 106 (1991), 139. 

[20] N. Rcshctikhin, A. Varchenko, Quasiclassical Asymptotics of Solutions to the KZ Equations [hep- 
th/9402126]. 

[21] H. Babujian, A. Lima-Santos and R.H. Poghossian, Int. J. Mod. Phys. A14 (1999) ,615 [solv- 
int/9804015]. 

[22] E. K. Sklyanin and T. Takebe, Phys. Lett. A219 (1996), 217 [q-alg/9601028]. 
[23] D. Bernard, Nucl. Phys. B309 (1988), 145. 

[24] A. Lima-Santos and W. Utiel, Nucl.Phys. B600 (2001), 512 [nlin.SI/0008002]. 

[25] V. Kurak and A. Lima-Santos, J. Phys. A : Math. Gen. 38 (2005), 333 [nlin.SI/0407041]. 

[26] V. Kurak and A. Lima-Santos, Nucl. .Phys. B701 (2004), 497 [nlin.SI/0406025]. 

[27] P. P. Kulish and N. Manojlovic, J. Math. Phys. 44 (2003), 676 [nlin.SI/0204037]. 

[28] E. K. Sklyanin, J. Phys. A: Math. Gen. 21 (1988), 2375. 

[29] I. V. Cherednik, Theor. Math. Phys. 61 (1984), 35. 

[30] K. Hikami, J. Phys. A: Math. Gen. 28 (1995), 4997. 

[31] A. Di Lorenzo, L. Amico, K. Hikami, A. Osterloh and G. Giaquinta, Nucl. Phys. B 644 (2002), 409 
[cond-mat/0206521]. 

[32] M. D. Gould, Y. -Z. Zhang and S. -Y. Zhao, Nucl.Phys. B630 (2002), 492 [nlin.SI/01 10038]. 

[33] M. D. Gould, Y. -Z. Zhang and S. -Y. Zhao, J. Phys. A: Math. Gen 35 (2002), 9381 [nlin.SI/0110038]. 



27 



[34] W. -L. Yang, Y. -Z. Zhang and M. D. Gould, Nucl. Phys. B698 (2004), 503 [hep-th/0401048]. 

[35] W. -L. Yang, Y. -Z. Zhang and R. Sasaki, Nucl. Phys. B729 (2005), 594 [hep-th/0507148]. 

[36] A. Lima-Santos and W. Utiel, Int. J. Mod. Phys. B20 (2006), 2175 [nlin.SI/0504060]. 

[37] L. Mezincescu and R. I. Nepomechie, J. Phys. A: Math. Gen. 24 (1991), L17. 

[38] H.-Q. Zhou, X.-Y. Ge, J. R. Links and M. D. Gould, Nucl. Phys. B546 (1999),779 [cond- 
mat/9809056]. 

[39] L. Mezincescu and R. L Nepomechie, Int. J. Mod. Phys. A7 (1992), 5657 [hep-th/9206047]. 
[40] V. V. Bazhanov and A. G. Schadrikov, Theor. Math. Phys. 73 (1989), 1302. 
[41] A. Lima-Santos, Nucl. Phys. B558 (1999), 637 [solv-int/9906003]. 
[42] G. -L. Li, K. -J. Shi and R. -H. Yue, Nucl. Phys. B670 (2003), 401. 

[43] S. Artz, L. Mezincescu and R. I. Nepomechie, J. Phys. A: Math. Gen. 28 (1995), 5131 [hep- 
th/9504085]. 

[44] R. E. Behrend, R A. Pearce and J. -B. Zuber, J. Phys. A: Math. Gen. 31 (1998), L763 [hep- 
th/9807142]. 

[45] V. Kurak and A. Lima-Santos, Nucl. .Phys. B699 (2004), 595 [nlin.SI/0406050]. 

[46] V. Kurak and A. Lima-Santos, J.Phys.A : Math.Gen. 38 (2005), 2359 [nlin.SI/0407006]. 



28 



A The Commutation Relations 



The equation (3.12) gives us the commutation relations for the matrix elements of the double-row mon- 
odromy matrix which play a fundamental role in the algebraic Bethe Ansatz. Here we present the 
commutations relations and their coefRcientes using a compact notation. Recall that the entries of the 
operator 

/ Viiu) Bi{u) B2{u) \ 
U{u) = Ci{u) V2{u) Bsiu) , (A.l) 
V C2{u) Caiu) Vsiu) J 

are given by (3.5). 

Substituting (2.15) and (A.l) into (3.12) we get 81 equations involving products of two matrix elements 
of U{u). These equations can be manipulated in ordem to put the product of pairs of operators in the 
normal ordered form. To do this we shall proceed in the following way. First we denote by E[i,j] = the 
component of the matrix equation (3.12) and collect them in blocks B[i, = 1, 5, j = i, ...,10— i, 
defined by 

B[t,j] = {F^J = E[i,jI /y = E[j, z], FFy = ^[10 - z, 10 - j], = ^[10 - j, 10 - i]} 

(A.2) 

Prom these blocks we can see that the pair {Fij,fij) as well as {FFij, ffij) can be solved simultaneously. 
We introduce the notation 

Di = Vi{u), di = Vi{v), Bi = B^{u), bi = B^{v), C\ = C^{u), Ci = Ci{v) (A.3) 

for the operators of the double-row monodromy matrix and 

Xi = Xi{u + v), Yi = yi{u + v), Xi = Xi{u-v), yi = yi{u-v), (A.4) 

for the Boltzmann weigths and 

{Z}i3 = = {^}ij('".^)> Fi = fi{u), fi = fi{v). (A.5) 

for the coefHcientes of the commutations relations, where 

{Z] = A,B,C,D,E,X and {z} = a,h,c,d,e,x. (A.6) 

Taking into account these simplifications, we will indicate the pair {Fij,fij) or {FFij, ffij) for which 
the corresponding normal ordered relations were obtained: 

Dibi = anbiDi + ai2Bidi + ai3Bid2 + auB2C\ + ai5S2C3 + ai6&2C'i 

Cirfi = AiidiCi + A12D1C1 + A13D2C1 + AuBiC2 + A15B3C2 + AiQbiC2 (A.7) 
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where the coefficients are 



{Fi7,fn) 



_ X1X2 . _ p ^5X2 , _ X5 

^11 - —ir^ Ai2 — -l<l— — , ^13--^, 

ynXe . _ ^7 4 _ 2:1 ^6 „^ 

^14 — — T^, ^15 — ^16 — -r^- 

X2 Xi Xi X2 Xi 



D162 = bnb2Di + bi2B2di + bi3B2d2 + bi4B2ds + bisBibi + biQBib^ 

C2di = BndiC2 + B12D1C2 + B13D2C2 + BUD3C2 + B15C1C1 + BieCsCi (A.9) 



where 



X1X3 ye Xq Xr yr X3 ye Xq Xr 

X3 Xi X3 Xi Xi X3 Xi xs Xi Xi 



Bi4 = 5,5 = -^^, ^i6 = -f. (A.10) 



Dibs = Xiik^Di + xi2biDi + xisBidi + Xi4Bid2 + Xi^Bids + XieB2Ci 

+X17B2C3 + xisb2Ci 
Cadi = XiidiCs + Xi2diCi + Xi^Dici + X14-D2C1 + X15D3C1 + X16-B1C2 

+XirBsC2 + XisbiC2 (A.ll) 

with the following coefficients 

^ 2:2X3 2/5 2/6 >5 ^ 

2:3 A2 X3 X2 X3 X2 

V p 2/6 ^4 p ^6 2/7 >6 2/6 , 
Ai3 = i^i — t2— -r; -r^^l2 

Xz X2 X2 X3 X2 X3 X2 

y p ^6 , 2/6 ^4 2/6 . 

X2 Xz X2 Xz X2 

Y ^6^2/7 2/6 f. 
-^15 = ^-T^j ^16 = T^^14 

A2 A2 
3^ 2/6 ^5 2/6 5^5 . ^ ys 2/6 5^5 . , . , „^ 

XZ X2 Xz X2 Xz Xi, A2 

Note that for each pair of equations the corresponding commutation relations are related by inter- 
changing 

Di ^ di, Bi ^ Ci, C'i ^ bi (A. 13) 
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• (-F24,/24) 

Cih = ciifeiCi + C1261C3 +C13B1C3 + C14-B3C3 + 01562(^2 +ci6C?i-Di +017^1 D2 

+ci8-Didi + cigDid2 + cnQD2di + Ciii£»2rf2 (A.14) 

with the following coefBcients 

X4 X5 Xq 2/5 Xq X-j X5 

Cll = C12 = — , Ci3 = — , Ci4 = --T^, Ci5 = — , 

Ai X2 Ai X2 Ai Ai Ai 

Y5 , IP X2 X5X2 t iV^ X2 ^ ^ Xr, 

C16 = -^ + Fi — — , C17 = — — , C18 = -h{ — ^+Fl^), 

Ai a;2Ai a;2 Ai a;2 Ai Ai 

ci9 = -( — l^+^i-^)' ciio = -/i — , ciii = - — . (A.15) 
a;2 Ai Ai Ai Ai 

The pairs (Fi6,/i6), (-Fis^/is), {FFie, f he), (^^i8,//i8) and F19 form a closed set of equations 

from which we have derived the following commutation relations 

-B261 = 61161^2 + 61262-81 + 61362-83, 

C1C2 = i^iiC2Cl + -E'l2CiC2 + -B13C3C2 

-B362 = 64162-83 + 64261-82 + 64363-82 

C2C3 = -B41C3C2 + E42C2CI + -B43C2C3 

-8162 = 62162-81 + 62262-83 + 62361 B2 + 62463-82 

C2C1 = E21C1C2 + E22C3C2 + E23C2C1 + i?24C2C'3 

-B263 = 63163^2 + 6326152 + 63362-Bl + 63462-B3 

C3C2 = -E31C2C3 + -E32C2C1 + i?33CiC2 + -E34C3C2 

.8262 = 62B2, C2C2 = C2C2 (A.16) 



where 



X2 X2 2/5 X2X6 

eil = — — , 612 = —, 613 = T^- (A.17) 

Xi A3 Xi Xi A3 

E41 = — -^42 = — -^43 = — (A.18) 

a:;i A3 xi A3 a;i 

2:2 X5{xj + x"^ - x^yr,) X2XQ 



xi Xl + Xgle ' 2^1 {xl - xzyr) X| -f XqYq 

~ xl- x^yr, XI + XeYe ^ x^ X| + ^6^6' ~ xl- x^yr, X| + X^Y^ ^"^^ 

X1X2 X2X3 a;5(xi + X2 - .xsj/s) X2YQ 



xj - X5y5 Xi + XeYe ' xi {xj - x^y^) X| + XqYq 

^ X2 X^jYQ ^ Xr, XqYq XlX^ X2 , , 

"""" ariXl+M' a^iXl+M - 2:5^5 X| + We ^ ^ ^ 



2 
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• (-F45,/45) 



D2bi = a2i?>i-D2 + 022-81 rfi + 023-81^2 + 024-83(^1 + a25-B3ci2 + a26-B2Ci 
+027-8203 + a28&2C'i + a^^b^Cz 

Cid2 = A2ld2Ci +^22-DlCi + yl23-D2Cl + A24-D1C3 +^25-D2C3 +^26-BlC2 

+A27-B3C2 + A2sbiC2 + A29b3C2 (A.21) 



with 



X4 X4 X4 Y5 Xq 

021 — -r; -^3 V TT^IZ -^14 

X2 ^2 X2 X2 X2 X2 X2 

f VbXi , y5 Yb. X4Y5 xe 
0.22 = fi — ^ - (i^i H ^)ai2 ^^11 ^12 

X2 X2 X2 X2 X2 X2 X2 

2/5 ^4 / 1;, , 2/5 ^5 s , Xq X^ 

023 = — -TP -rr/On-, a24 = ji-rp ^11, 025 = -r;^ 

X2 X2 X2 X2 X2 X2 X2 

. p , yS >5 ^ X4Y5 Xq 

026 = -[Fi-\ ^)ai4 ^^16 -'^18 

X2 X2 X2 X2 X2 

2/5 ^5 / , 2/5 ^5 X 2/6 / r , 2/5 >5 . X4 Y5 Xq 

027 = -TP {Fl-\ -r^)Ol5 , 028 = {Fi-\ T^jOia T^^14 A16 

X2 X2 X2 X2 X2 X2 X2 X2 X2 X2 

X4X5 X4Y5 Xe 

029 — -Tz ^^15 (A.ZZ) 

X2 X2 X2 X2 X2 



(-^15) /15) 



-8161 = 60161-81+60262-02 + 60362-01 + 604-82(^1 + 605-82^2 

CiCi = i;oiClCi+-E;o2rf2C2 + -Bo3rflC2 + -Bo4-DiC2 + -Bo5-D2C2 (A.23) 



with 



a;3a:;4 - xaye X3X4 - x^ye Xe 

eoi = , eo2 = — 

X1X3 X1X3 X2 

xsVe - xeVr X3 X3X4 - xeye Xq 
eo3 = -T^ + Fi -TT- 

X1X3 X2 XiXs X2 

Xe X3 Xe Xq 
eo4 = - — ^-bfi — , eo5 = ^ (A.24) 

X3 A.2 yi.2 -^2 

These commutaion relation play a special role in the construction of the n-particle states, as men- 
tioned above. 

• {FF2e,ff26) 

-B163 = ^1163.81 + ^1261-81 + rfi362-Dl + ^1462-02 + 615-82^1 + (^16-82^2 + (^17-82^3 
C3C1 = -DiiCiC3 + -Di2CiCi+-Di3rfiC2 + £'l4rf2C2 + £115-0162 + £'l6-D2C2 + -Di7-D3C2 

(A.25) 
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where 



and 



where 



Di3 

Die 



Y^Bie-X^' ,T;2(nBi6 -^4) 

{yr,Y„ + .r2>5^i^j£(i:i - /i/y^-y2 + ■'■2Yr,Bn 

X2{Y^Bie-Xi) 
{y^Ye + a:2y5-Bi5)-go2 - t/s^a 

xaCFsBie -X4) 
(ysVe + X2Y5Bi5)Eo4 + ^1^5X2 + F2X2X5 + X2Y5B12 

x2{Y^Bie-Xi) 
{yhYe + X2Y^Bi^)Ea^ + 775^2 + j^3.x2^5 + 3:2^5^13 
X2{Y^Bie-X4) 



^ir = -^i^ (A.26) 



-B3&I = ^2161-83 + rf22?>l-Bi + rf23?>2-Dl + d24?>2-D2 + rf25''2-D3 + ^26-82^1 + (^27-82^2 
C1C3 = £'2lC3Cl + -D22C1C1 + D23d\C2 + D2id2C2 + D25d3C2 + -D26-D1C2 + D27D2C2 

(A.27) 



^ ^4 ^ ^5 R , ?/5 1^6 Ye Y. 

«21 = -T^ + T^-D16, «22 = — T^eoi + — £>15 

A3 A3 a;2 A3 a:;2 A3 A3 

c^23 = F,yi^ + F2^-'-^^eo3 + ^B,2 

X2 A3 A3 X2 A3 A3 

ys ^2 ^ p ^5 ^6 ^ n o , ^5 ^ n p 

«24 = -TT + -f'STF -Tr^02 + -r^-Dl3, «25 = + ^-Ol4 

X2 A3 A3 a;2 A3 A3 A3 A3 

, r ?y5 ^2 ?y5 Ye Y^ 
"26 = -Ji — -r^ T^eo4 + — i>ii, 

2:2 A3 X2 A3 A3 

rf,, = _^^_^^eo5 (A.28) 



(-F28) /28) 



X2 X3 X2 Xs 

-D2&2 = 62ii>2-D2 + 622^2^1 + 623-62^2 + 624-82^3 + 625-B161 + 626-8163 

+627-B361 + 628-B363 
(^2^2 = B2ld2C2 + B22D1C2 + B23D2C2 + B24D3C2 + B25C1C1 + B26C3C1 

+B27C1C3 + B28C3C3 (A.29) 
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where 

-D21 = J- H -r^"27 + ^eo5, -D22 = — /1-D12 H -r^«23 + ^603 

a;2 A2 A2 a;2 A2 A2 

-D23 = — /I-DI3 H "r^"24 + -r^eo2, -t>24 = — /l-t'14 H "^"25 

a;2 A2 A2 X2 A2 

D ^ D , ^5 -'sTa Ye „ ^ d , ^5 , 

-D25 = — 7li>15 H -r^"22 + ^eoi, -t>26 = -/1-D16 H 1^"21 

X2 A2 A2 a;2 A2 

^27 = -^f^, 528 = ^ (A.30) 

a;2 A2 A2 

(-F-P35, //as) 

C3&1 = C2161C1 + C2261C3 + C23-B1C3 + C24-B3C3 + C2562C2 + 026^1 r>i 

+027^1-02 + 028^1 -D3 + C29-Dl(ii + C210-Dld2 + 0211-02^1 + 0212-02^2 
+C213-D3dl + 0214-03^2 
C163 = C2161C1 + C22&3C'l + + C24-B3C3 + C2hi>2C2 + C26diDi 

+C27d2Di + C2»dsDi + C29Didi + C2wD2di + C2iiDid2 

+C213Did3 + C2uD2d3 (A.31) 

where 

2/6^5/ 1 ^ 2:4 2/6 yr , 2/6 

C21 = — -r^(cii-l), C22 = 1 -rrCi2, C23 = 1 -rrCi3 

X3 A2 X3 X3 X2 X3 X3 X2 

ye Xr, yQ Y5 ye Xi ye Y5 

C24 = T^Ci4, C25 = — H -Trci5 

X3 X2 X3 X2 X3 X2 X3 X2 

X4 Fe „ X6 X3 ye Y7 ye Y5 

C26 = -fl -^2 H t;^C16 

X3 A2 a;3 A2 a;3 A2 0:3 A2 

rre X3 2:4 Ye , ye >5 are X3 

C27 = --^a — -rr ^ -^r ^ ^^17, C28 = 

X3 X2 X3 X2 X3 X2 X3 X2 

f yr Ye ye Xi Xe ye Y5 

C29 = fl — r^-flFi — — + fiF2—-\ ^Ci8 

X3 X2 X3 X2 X2 X3 X2 

2/7 Ye ye X4 Xe ye Y5 
C210 = — -r; — -^ + ^2—^ -r^Cig, 

X3 X2 X3 X2 X2 X3 X2 

2/6 ^4 , r Xe ye Y5 

C211 = -/i — -r^ + /1^3-r^ H ^ciio 

X3 X2 X2 X3 X2 

ye X4 Xe ye Y^ Xq Xq 

C212 = — + i'3T^H T^Clll, C213 = h^, C2IA = ^ (A.32) 

2:3 A2 A2 2:3 A2 A2 A2 

(Fi^2a, //23) 

-D361 = a3i&i-D3 + 032-61^1 + 033-81(^2 + 034-83(^1 + 03553(^2 + a36-B2Ci 
+a37B2C3 + a2sh2Ci + 02962(73 

Cxd3 = A3id3Ci + A32DXC1 + ^3a-D2Cl + A3iDxC3 + A3r,D2C3 + ^36^102 

+A37B3C2 + A2861C2 + A2963C2 (A.33) 
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where 



X2 X2 X2 Yd 

0-31 = — -r; 

2^3 -'^3 2:3 A3 

132 = -/I -TT-^ -TT^ll - {Ql)0'22 - (Q2)ai2 -^T^^'^ 

X3 A3 X3 A3 X3 A3 

033 = - — ^ - ((32)ai3 - (Qi)a23 
X3 A3 

„ ye X2 . , x2Y(i ye X2 , 

0134 = Ji — -TT (^1)024 -rrXii, 035 = — -7? [Qi)a25 

X3 A3 X3 A3 0:3 A3 

0136 = —-^Aie - {Qi)a26 - {Q2)ai4 - —-^Xis 
X3X3 X3X3 

037 = - — ^ - (<32)ai5 - (Ql)a27 

X3 A3 

038 = — ^^14 - (Q2)ai6 - (<5i)a28 - — 

X3 A3 X3 A3 

039 = — T^H -r^^i5 - (Qi)a29 -^Xir (A.34) 

2:3 A3 X3 X3 X3 A3 



{FFi3,ff 13) 



-D3&2 = 631-82(^3 + &32-B2dl + 633^2^2 + &34&2-D3 + 635S161 + 636-B163 

+637-6361 + 638-B3&3 
(72^3 = i331-D3C2 + i332£>iC2 + -B33-D2C2 + 534^3(^2 + S35C'lCl + S36C3C1 

+S37C1C3 + ^38^303 (A.35) 



where 



yr Xi 
X3 X3 

—f2—-rr H — --^Di3 H — --J-{Bii + Dis^ie + £'03-815) 
2:3 -'^3 X3 A3 0:3 A3 

-((32)612 - (Ql)622 



631 = - — T;i-(g2)6i4-(gi)624 



632 



633 = -/3 — ^ + —^Z)l4 + — ^(^4-816+ £02-615) -(Q2)6l3-((3l)623 

2:3 A3 Xs A3 2:3 A3 

634 = ^^ + ^^iB,, + D,rB,e) + ^^Dn 

X3 X3 X3 X3 X3 A3 

635 = —^Di2 + — ^(£'l2Bl6 + -SoiSls) - (Q2)6l5 - (Ql)625 

X3 -^3 X3 A3 

636 = (02)616 -(Ql)626 

X3 X3 

637 = -^I)ll + -^i)ll-Bi6-(gi)627 

628 = --^-(Qi)628 (A.36) 
X3 A3 
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Here we have used the notation 

Q,=F3-^f, Q. = F.-F,^f + (A.37) 

X3 A3 2:3 A3 2:3 A3 

B The H amplitudes 

In this appendix we summarize the expressions for the i?-functions presented in (3.18) and (3.25) 

Hn{Up,Uq) = ai4{u,Up) (ci6(Up,Wg) + Cl8{Up,Uq)) +ai5{u,Up) {C26{Up,Uq) + C29{Up,Uq)) 

+6i2(u, Up)Gdi {up, Uq) + uj{up, u}an{u, Up)ai2{u, Uq)Gd^ {u, Up) 

Hi2{Up,Uq) = aii{u,Up) {ci7{Up,Uq) + Cuo{Up,Uq)) + au{u,Up) (C27(Wp, Wg) + C21l{Up,Uq)) 

+6i3(w, Up)Gdi (up, Uq) + u>{up, u)an (u, Up)ai2{u, Uq)Gd^ (u, Up) 

Hl3{Up,Uq) = ai4{u,Up)ci^iUp,Ug) + ai5{u,Up)c21(:)iUp,Uq) + bl2{u,Up)Gd2iUp,Uq) 

+U!{up, u)aii{u, Up)ai3{u, Uq)Gdi i'^i 

Hu{Up,Uq) = ai4iu,Up)ciii{Up,Uq) + ai5iu,Up)c212{Up,Uq) + biz{u,Up)Gd^{Up,Uq) 

+uj{up, w)oii (u, Up)ai3{u, Uq)Gd2 (u, Up) (B. 1) 

and 

Hal{Up, Uq) = aje{u, Up) (ci6(Mp, Uq) + CisiUp, Uq)) + aj7{u, Up) {C2&{up, Uq) + C29{Up, Uq)) 

+bj2{u, Up)Gdi {up, Uq) + u;{up, u)aji{u, Up)aj2{u, Uq)Gdi {u, Up) 
+aji {u, Up)aj4,{u, Uq)di3{up, u) 

Ha2{Up,Uq) = aj6{u,Up) {cn{Up,Uq) + Ciio(%,Ug)) + aj-j{u,Up) {C27{Up,Uq) + C21l{Up,Uq)) 
+bj3{u, Up)Gdi {Up, Uq) + U){Up, U)aji{u, Up)aj2{u, Uq)Gd2 {u, Up) 

+aji{u, Up)aj4,{u, Uq)di4{up, u) 

HasiUp^Uq) = aje{u,Up)cig{Up,Uq) + ajjiu,Up)c21oiUp,Uq) + bj2{u,Up)Gd2iUp,Uq) 

+uj{up, u)aji (u, Up)aj3{u, Uq)Gd^ {u, Up) + aji (u, Up)aj5 (u, Uq)di3{up, u) 

Hai{Up,Uq) = ajeiu,Up)cilliUp,Uq) + aj7{u,Up)c212{Up,Uq) + bj3iu,Up)Gd2iUp,Uq) 

+Lo{up, u)aji {u, Up)aj3 {u, Uq)Gd2 (w, Up) + Oji (w, Up)aj^ {u, Uq)dii{up, u) (B.2) 

for a = 2, 3. 
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These function satisfy the properties 

Hal{Up,Uq) = ij{Up,Uq)Hal{Uq,Up) 

Ha2{Up,Ug) = Uj{Up,Uq)Has{Uq,Up) 

Ha3{Up,Uq) = ij{Up,Uq)Ha2iUq,Up) 

Hai{Up,Uq) = Ul{Up,Uq)Ha4{Uq,Up) (B.3) 

and in the quasi-classical limit of (5.5) their expansion have to be considered up to rj'^ order: 

Haj {up, Uq) = h'f] {up, Uq) + h')^^ {up, Uq)7] + ^h^^^ {up, Uq)rf + ^k^^] {up, Ug)?]^ + o{r]'^) 

a = 1,2,3, j = l,2,3,4 (B.4) 

where h'f^{up,Uq) = h^^]{up,Uq) = 0. 
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